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We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and 7— ray wavelengths 
of the BL Lacertae object AO 02 35+164 in 2008. H ere we reproduce the results of the analysis of these 
observations, already presented in lAgudo et al.l [201 lbjl ■ We analyze the timing of multi- waveband variations 
in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7 mm 
with ~ 0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at 
high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks 
in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a super luminal 
jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to 
conclude that the outburst occurred in the jet both in the quasi-stationary core and in the superluminal knot, 
both at > 12 parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence 
of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes 
through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. 
The multi- wavelength light curves vary together on long time-scales (months / years) , but the correspondence is 
poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the 
outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role 
in modulating the synchrotron and inverse Compton fluxes. 



1. Introduction 



Our understanding of the processes leading to 
the generation of 7-ray emission from blazars, the 
most extreme active galactic nuclei, depends on 
where those 7-rays originate. This is currently 
the su bj ect of considerable d e bate [e .g.. lAgudo et al. 
201 lal iMarscher fc Jorstadl l2010l iTavecchio et al.l 
2010] . Two main locations of the site of 7-ray emis- 
sion in blazars have been proposed. The first is close 
(< 0.1 — lpc) to the supermassive black hole (BH). 
The second one, a region much farther (>> lpc) 
from the BH, is beyond the "core" where the jet 
starts to be visible at millimeter wavel engths with- 



3C 4 54.3, respectively. 



VLB! Supporting this second scenario, lAgudo et al.l 
[2011a I unambiguously locate the region of 7-ray flares 
> 14 pc from the BH in the jet of OJ287 through 
correlation of millimeter-wave with 7-ray light curves 
and direct ultrahigh- resolution 7 mm i maging with the 
VLBA . Si milar results are obta ined by IMarscher et all 
[2010] and lJorstad et al.l [2010] for PKS 1510-089 and 



In lAgudo et al.l [2011bl ] we investigate the location 
and properties of a radio to 7-ray outburst in the 
BL Lacertae object AO 0235+164 (0235+164 here- 
after, z — 0.94), which results are reproduced here. 



2. Observations 

Our photo-polarimetric monitoring observations of 
0235+164 (Figs. [T][3]) include (1) 7 mm images with the 
VLBA from the Boston University monthly blazar- 
monitoring program, (2) 3 mm observations with the 
IRAM 30m Telescope, and (3) optical measurements 
from the following telescopes: Calar Alto (2.2m 
Telescope, observations under the MAPCAT pro- 
gram), Steward Observatory (2.3 and 1.54m Tele- 
scopes), Lowell Observatory (1.83m Perkins Tele- 
scope), San Pedro Martir Observatory (0.84m Tele- 
scope), Crimean Astrophysical Observatory (0.7m 
Telescope), and St. Petersburg State University (0.4m 
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Figure 1: Sequence of 7 mm VLBA images of 0235+164 convolved with a FWHM = 0.15 mas circular Gaussian beam. 
Contour levels represent total intensity (levels in factors of 2 from 0.4 to 51.2% plus 90.0% of peak= 4.93 Jy/beam), 
color scale indicates polarized intensity, and superimposed sticks show the orientation of \- Reproduced from 
lAgudo et all [2011bl ]. 



Telescope). Our total flux light curves (Fig. [2]) include 
data from the Fermi-LAT 7-ray (0.1-200 GeV) and 
Swift-XKT X-ray (2.4-10 keV) observatories, avail- 
able from the archives of these missions, and RXTE 
at 2.4-10 keV. Optical i?-band fluxes come from the 
Tuorla Blazar Monitoring Program, the Yale Univer- 
sity SMARTS program, and Maria Mitchell and Abas- 
tumani Observatories. Longer wavelength light-curves 
were acquired from the Submillimcter Array (SMA) at 
850 /im and 1 mm, the IRAM 30m Telescope at 1 mm, 
the Metsahovi 14m Telescope at 8 mm, and both the 
Owens Valley Radio Observatory (OVRO) 40m Tele- 
scope Fermi Blazar Monitoring Program and Univer- 
sity of Michigan Radio Astronomy Observatory (UM- 
RAO) 26m Telescope at 2 cm. 

We followed data reduction procedures de s cribed 
in previous studies: VLBA 
optical polarim etric data: 
IRAM data: lAgudo et all 
iGurwell et ail 120071 : Metsahovi 



Jorstad et al.l 
iJorstad et al.| 

~3L 



2006, 



Marscher et al 




201L 
20111 



2010f 



2001 
201C 



20K 



SMA: 



Terasranta et alj 



19981: OVRO: iRichards et all l201ll; U MRAO 



Aller et al.Hl985fl: Sw ift: I.Torstad et al.ll2010»: RXTE: 



an d Fer mi-LAT : 
Marscher et~aH [2010j . 



See 



Agud o et al 
Agudo et al 



for further details about the X-ray and 7-ray 



ever detected in this object) propagating from the core 
at < /3 app >= (12.6 ± 1.2) c. 

The ejection of Qs was coincident with the core near 
the start of an extreme mm outburst (08 mm in Fig. [2])). 
Figure [2] shows that radio and mm outbursts in 2008 
(08 ra d and 08 mm ) contain contributions from both the 
core and Qs. Their contemporaneous co-evolution 
suggests that the disturbance responsible for the ejec- 
tion of Qs extended from the location of the core to 
Qs in the frame of the observer, w hich could have re- 



sulted from li ght-travel delays [e.g., Agudo et alj200ll 
~l ll997l iMimicaet al j|2009j . Theraritv of 



I Gomez et al 



data reduction. 



the 08 ra d, 08 mm , and Qs events strongly implies that 
they are physically related. 

The je t half-opening-angl e of 0235+164 [a- mt /2 < 
1°25, see lAgudo et all2011bj [ and the average FWHM 
of the core measured from our 31 VLBA observ- 
ing epochs in [2007,2010] ((FWHM coro ) = (0.054 ± 
0.018) mas), constrain the 7mm core to be at d C ore = 
1.8(FWHM coro )/tana in t > 12pc from the vertex of 
the jet cone. 



4. Contemporaneous Flares from 7-ray 
to Radio Wavelengths 



3. Major Millimeter Flare in 2008 Related 
to a New Superluminal Knot 

We model the brightness distribution of the source 
at 7 mm with a small number of circular Gaussian 
components (Fig. |TJ). Our model fits include a bright 
superluminal feature (Qs, the brightest jet feature 



Figure [3] reveals that the 08 ra d and 08 mm flares 
were accompanied by sharp optical, X-ray, and 7-ray 
counterparts (08 opt , 08x, and 08 7 flares, respectively). 
Our formal light-cur ve correlation analys is (Fig. 2]) 
-performed following lAgudo et al.l [2011a — confirms 
the association of 7-ray variability with that at 2 cm, 
8 mm, 1 mm, and optical wavelengths at > 99.7 % con- 
fidence. The flux evolution of the VLBI core is also 
correlated with the 7-ray light curve at > 99.7 % con- 
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Figure 2: Left: Light curves of 0235+164 from 7-ray to millimeter wavelengths. X-ray photon index evolution from the 
Swift-XKF data is also plotted. Vertical dotted lines mark the three most prominent 08 op t optical peaks. The yellow 
area represents the time of ejection of feature Qs within its uncertain t y. RJD = Julian Date — 2400000.0. Right: Same 
as left panel for RJDe [54500, 55000]. Reproduced from lAgudo et al.l j2011bl |. 



Science. Moreover, the evolution of the degree of op- 
tical linear polarization (p pt) and X-ray light curve 
are also correlated with the optical i?-band, 1 mm, 
and 2 cm light curves at > 99.7% confidence (Fig. [5]), 
further indicating that the extreme flaring activity re- 
vealed by our light curves is physically related at all 
wavebands from radio to 7-rays. 

There is, however, no common pattern to the dis- 
crete correlation function (DCF) at all spectral ranges. 
This implies that, although there is correlation on long 
time-scales (years), on short time-scales (< 2 months) 
the variability pattern does not correspond as closely. 
This is the result of the intrinsic variability pattern 
rather than the irregular time sampling at some spec- 
tral ranges. 

5. Correlated Variability of Linear 
Polarization 

Figure [3] reveals extremely high, variable optical po- 
larization, p opt > 30%, during the sharp 08 op t opti- 
cal peaks. Whereas the integrated millimeter-wave 



degree of linear polarization (p mm ) and that of the 
7 mm core remain at moderate levels < 5 %, the po- 
larization of Qs (pmm.Qs) peaks at the high value of 
~ 16% close to the time of the second sharp opti- 
cal sub-flare. The coincidence of this sharp maximum 
of Pmm.Qs in the brightest superluminal feature ever 
detected in 0235+164 with the (1) high optical flux 
and polarization, (2) flares across the other spectral 
regimes, and (3) flare in the 7-mm VLBI core, implies 
that the ejection and propagation of Qs in 0235+164's 
jet is physically tied to the total flux and polarization 
variations from radio to 7-rays. 

On long time-scales (years), the linear polarization 
angle at both optical (Xopt) and millimeter (xmm and 
Xmm) wavelengths varies wildly, without a preferred 
orientation or systematic common trend. However, 
during flare 08 op t, Xopt maintains a stable orientation 
at (100 + 20)°, whereas xSm i s roughly perpendicular 
to this (~ 0°), as expected for a plane-perpendicular 
shock wave propagating to the south towards Qs. Ow- 
ing to the large peak value of Qs, p^rnQs ~ 16%, 
one cannot explain the orthogonal optical-millimeter 
polarizations by opacity effects. Instead, we propose 
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Figure 3: Left: Long term optical and millimeter- wave linear polarization evolution of 0235+164 in the 

RJD= [54000,55600] range. Right: Same as left panel for RJDe [54530,54850]. Reproduced from lAgudo et al.1 [2011bl |. 



that the optical polarization mainly arises in a con- 
ical shock associated with the 7-mm core, while the 
millimeter-wave polarization results from a propagat- 
ing shock front associated with Qs. We surmise that 
the moving shock also emits polarized optical radia- 
tion, since the optical polarization drops precipitously 
when the orthogonal polarization of Qs peaks (Fig. [SJ- 
right). 



6. Low Probability of Chance 
Coincidences 



7. Observational Evidences 

The coincidence of the ejection and propagation 
of Qs -by far the brightest non-core feature ever re- 
ported in 0235+164- with the prominent 7-ray to ra- 
dio outbursts and the extremely high values of p pt 
and Pmm.Qs provides convincing evidence that all these 
events are physically connected. This is supported by 
probability arguments and by our formal DCF analy- 
sis, which unambiguously confirms the relation of the 
7-ray outburst in late 2008 with those in the opti- 
cal, millimeter-wave (including the 7-mm VLBI core) 
and radio regimes. We locate the millimeter core at 
encore <^ 12 pc from the vertex of the jet. 



The relationship among the 7-ray, optical, and 
radio-millimeter flares is supported by probability ar- 
guments. If the flares occur randomly, the probability 
of observing, at any time, a 7-ray outburst like the one 
reported here (i.e., with flux > 10 -6 phot cm~ 2 s _1 
and duration ~ 70 days) is p 7 = 0.08. For optical 
and radio-millimeter wavelengths, this probability is 
Popt = 0.04 and p mm — 0.15, respectively. Thus, if 
the flares at different wavelengths were random and 
independent of each other, the probability of observ- 
ing a 7-ray, optical, and radio-millimeter flare at any 
given time is p 7j0 pt,mm = 5 x 10~ 4 . This counters 
the null hypothesis of random coincidence at 99.95% 
confidence. 



8. Conclusions 

We identify the 7-mm core as the first re-collimation 
shock near the en d of the jet's acceleration and colli- 
mation zone [ACZ ljorstad et al~ll2007l iMarscher et al.l 
l2008l l2010j . Superluminal feature Qs is consistent 
with a moving shock oriented transverse to the jet 
axis, given the extremely high p m m,Qs, with xSm P ar_ 
allcl to the direction of propagation of Qs. The flux 
evolution of the core appears closely tied to that of Qs, 
and its light curve is correlated at high confidence with 
those at 7-ray, optical, and millimeter wavelengths. 
This suggests that Qs is the head of an extended dis- 
turbance, perhaps containing a front-back structure 
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Figure 4: Grid of DCF of labeled light-curve pairs during the maximum time period RJD = [52200, 55600]. Top row of 
panels show DFC with 7-ray light-curve, whereas bottom row show DCF with .R-band light-curve. Grey dotted curves 
at positive (negative) DCF values symbolize 99.7% confidence limits for correlation against the null hypothesis of 
stochastic variability . Green dashed lines at and —50 day time-lags are drawn for reference. Reproduced from 
lAgudo et all [2011b| . 
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Figure 5: Same as F ig. [4]but for the 1mm light-curve (top) and the 2 cm light-curve (bottom). Reproduced from 
lAgudo et all [2011bl ]. 



stretched by light-travel de lays in the observer's frame 
[see, e.g.. lAlov et~aT1l2003{ . 

Under this scenario, the radio/millimeter- wave and 
optical (and perhaps X-ray) synchrotron flares start 
when the front region crosses the conical shock at 
the core, where the jet is at least partially optically 
thin (Fig. [5]). This interaction accelerates electrons, 
and produces inverse Compton 7-ray emission from 
the up-scattering of IR-optical photons. When the 
back region of the moving perturbation encounters 
the core, their interaction again produces efficient par- 
ticle acceleration, which is seen as a sudden optical 
and radio/millimeter synchrotron emission enhance- 
ment. The subsequent optical variability is produced 
by the passage of the remaining shocked turbulent 



plasma in the back structure through the core. Dur- 
ing the different optical sub-flares, the integrated ra- 
dio/millimeter synchrotron flux keeps rising. This ra- 
dio/millimeter outburst is more prolonged owing to 
the longer synchrotron cooling-time of electrons radi- 
ating at these wavelengths and the lower speed of the 
back structure. Indeed, Qs does not reach its maxi- 
mum radio / millimeter- wave flux until traveling a pro- 
jected distance of ~ 0.13 mas from the core. When 
the entire front-back structure passes across the core, 
the synchrotron emission declines rapidly at optical 
(and, if relevant, X-ray) frequencies, as does the 7- 
ray emission. The decay of the radio/millimeter- wave 
emission is more gradual (see above). 
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Figure 6: Scheme of proposed model for the multi-wavelength flaring behavior of AO 0235+f64. Radio-loud AGN 
sketch adapted from iMarscherl [ 2006f ] . 
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